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1 
2 
In the familiar Fibonacci sequence, 3 
the units position has been marked off. This 
sub-sequence repeats on a cycle of 60, as can be 5 
readily verified in a few minutes. 
8 
Let us generalize this single-digit sequence 
as follows: 13 
2} 1 
i aed @ starting values 3 4 
3 8/9 
6 1 44 
Be 20 2 3)3 
9 7 
5 T 
2 4 
1 
5 
6 
starting with two columns, both initialized to one, and i 
with one as the next term of the first column. en 7 
each successive term is formed by adding the preceding 
two terms, alternating columns as shown, Only one 8 
digit is retained in each addition (i.e, all arithmetic 
is modulo 10). 5 
This pattern, too, must repeat, and it does, on 3 
a cycle of 217 terms (that is, 217 rows of the 2-column 
@ pattern). 8 
1 


Continued on page 3 
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Write for 


’] 

Log 40 1 
In 40 3 
Vio 6 
Vio 3 
Vio 1 

'ViiO 

et 

ho 


SOMEWHERE IN HERE::,; 


there is an article 
about computers 
that may be 

helpful to you. 


These are the 203 
business and 
computer journals 
we read 
every month... 
and digest the best 
of the articles 


for 
YOU 


information: 


If you need to 
keep up with 

the computer field 
but lack the 

time to research 
and read, 


DATA 
PROCESSING 
DIGEST 


is your answer. 
Monthly, averaging 
50 items, including 
book reviews; index, 
calendar, complete 
references to 
original articles. 
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Extending this scheme to three columns, we have: 


starting ( Se 


values il 2 3 
a 5 
o 64 
CONTEST 
6 6 0) Ol LO) 
5 1 
OF 5 
6 7 
2 po5 
Pee ale” © 2 


and the 3-column pattern is found to repeat on a cycle 
of 520 terms. The cycle lengths for the 4- and 5- 
column patterns are readily obtained, and we have the 
following table: 


Number of Cycle 
columns, K length 


60 
217 


520 
42 
196812 


indicating that we have an irregular (indeed, weird) 
function. 


Our 9th contest, then, will award our usual $25 
prize to the person who extends table T the furthest. 
The computer program used must be furnished. All entries 
must be received by POPULAR COMPUTING, Box 272, Calabasas, 
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K-COLUMN FIBONACCI 


California 91302, by September 30, 1976. fel 
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Using New Tools 


_ ART-.: COMPUTING 10—# 


When Numerical Methods are devised to implement 
the techniques of Numerical Analysis, troubles arise 
from several sources: 


1. In a computer, most numbers do not exist. In 
a binary machine, the numbers .1, pi, and the square root 
of 5 do not exist, for example. If 8-digit scientific 
notation is used, there are just 10° numbers between zero 
and one that can be expressed exactly; all others are 
approximations. Multiple precision does not solve this 
difficulty, although it may relieve it. 


2. The numbers that do exist in a machine (again 
in scientific notation) are not uniformly dense. There 


are also just 10“ numbers between 10° and 107, which means 
that in that range the numbers are spread out thinner. 


3. The notion that A+B = BtA doesn't always hold. 
Consider the addition (in a 3-digit system) of: 


100. 
oa 


a million of these 


Hh ee 


. 
. 
. 
. 


If the addition is done from the top down, the result is 
100. If it is done from the bottom up, the result is 
an overflow. 


4. We can lose significance at any time, with no 
warning. The addition of: 


34567891 E06 
-34561234 E06 


(both correct to 8 significant digits) will be 

66570000 E02. 
The loss of significance is due entirely to the fact 
that the numbers are close to each other, but the point 
is that that loss takes place unseen. 


5. Due to the finite nature of all numbers in a 
computer, problems that are mathematically sound may be 


computationally unstable. The value of the determinant: 
1 2 3 
4 5 6 


10 20 30 


is clearly zero, since the first and third rows are 
proportional. Mathematically, the determinant: 


1 2 3 
4 5 6 
10 20 29.999999 
is non-singular, but in a computer it is "almost singular,' 


which may cause serious trouble when the situation is 
better disguised than it is here. 


6. Rounding errors do not always balance out. 
See, for example, the calculations on the perimeters of 
circumscribed polygons given in Richard Hamming's article 
"Archimedes and the Value of Pi" in our issue No. 12. 
Or, see the evaluations of the Taylor series for sine 
in Numerical Methods and Fortran Programming (McCracken 
and Dorn, page 89) in which the value for sine of 2190° 
is given by the obvious Fortran program as 25902480. 
Further interesting examples are found in the article 
"Some Dangers of Machine Calculations" by Leon Winslow 
in the Journal of Recreational Mathematics, Vol. 8, No. 2, 
page 83. 


Despite all this, most numerical processes do work 
and therein lies the real difficulty. When one applies 
a numerical process to a problem with a known (analytic) 
result, and that result is obtained (involving thousands 
of individual calculations), the troubles and pitfalls 
listed above seem to be remote; that is, they seem to 
apply only to unusual cases which can't happen to me. 


» 
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The Error Amplification problem (in Issue No. 32) 
was designed specifically to reveal such troubles--but 
it didn't. The expression: 


eet O.. A 1S 15 287 Pie Ong 99 


9 2 4 6 8 10 12 14 16 a aa 98 


was to be evaluated by many distinct methods (and, in 
particular, by taking the various powers and roots in 
sequence), with the thought that each different method 


would yield a distinct result. The result was somewhat 
the opposite; all the approaches yielded essentially the 
same result, within the limits of precision used. For 


example, the "true" result may be taken as 248.81398019578 
(done with 100-digit arithmetic, working first on the 

long exponent and then calculating the power by logarithms). 
Carrying out the calculation sequentially on a pocket 
calculator, using 12-digit arithmetic, the result is 

248 .8139762. The result has, if anything, reinforced 

our intuitive belief that all is well with the world. 
Providence again seems to be at work guiding fools, drunks, 
and numerical workers. 


Let us explore the theories involved here by doing 
some numerical integration by the method probably most 
widely used; namely, the formula of Simpson. 


Simpson's Rule provides a widely used method for 
numeric evaluation of integrals. The rule is: 


I= [roe ~ 2 }¥0 Hes Wis ik ACY Saq) + el 
‘a 
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where the interval from a to b is divided into an even 

number of sub-intervals of width h, and the y values in 
the formula are the values of the function f(x). The 

theory assures us of an exact value for polynomials of 

degree three or less. 


As with any new tool, we try it out first on known 
cases. Thus, for the curve 


y = x3 - llx® + 4x + 60 


the area between -1 and +3 should be obtainable precisely 
by Simpson's Rule using any even number of sub-intervals, 
and the result should agree with the analytic solution: 


xt 7/4 -(11x3)/3 + 2x? + 60x 
evaluated at +3 and -1; the result is 173 1/3 square units. 


We can apply the Rule with just two intervals (N = 2), and 
calculate 


_ 3- (-1) [ f 
I= Rr ae jay + 0 + 4(54)| 


where h = 4/2; 44 is the value of the function at x = -1; 
O is the value of the function at +3; and 54 is the value 
of the function at the midpoint of the range (where x = +1). 


The results agree; our new tool tests out exactly 
for a known case. We can proceed to try it out in other 
situations to see how powerful it is. For example, we 
can construct a 4th degree curve with roots at -1, 1, 9, 
and 10: 


y = (x+1)(x-1)(x-9)(x-10) [ 
y = x# ~ 19x3 + 89x" + 19x - 90 


PCc4O-8 


and find the area of the large arch (between 1 and 9) 

analytically to be 2286.9333. Then we can apply 

Simpson's Rule to the integral, using 2 intervals, then 

4 intervals, 6 intervals, 8 intervals, and so on, to find $ 
at what point the process gives a reasonable approximation 


to the desired area. The results are as follows: 

Number of Value by 
intervals Simpson 

2 2560 .00000000 

4 2304 .00000000 

6 2290, 30452681 

8 2288 .00000000 

10 2287 . 37024021 

20 2286 .9606 3983 

30 2286 .93872714 

50 2286 93403232 

70 2286 .93351483 

100 2286 .93337786 

120 2286 .93335497 


Again, we have found that our new tool works as it 
should. We will try it once more on a non-polynomial. 
We will calculate the area of a quarter circle of radius 
one, Which is given by: 


1 


[i - x? ax 
° 


whose value is 7/4 = .785398163397... 


Number of Value by 
intervals Simpson 

2 - 74401694 

ut -77089879 

8 . 78029729 

16 -78359942 

32 . 784.76 305 

64. -78517377 

128 . 78531885 

256 78537013 

512 - 78538825 
1024 -785 39466 . 

2048 -78539692 

4096 -78539773 


8192 - 78539801 


Having done all this, we should be convinced that 
our new tool provides a workable method for numerical 
integration. For an unknown integral, the only sticky 
problem is the number of sub-intervals to use to insure 
the level of precision we seek. Let us now apply the 
tool to one more integral: 


whose value is readily found to be exactly 8. 


The 
following are some preliminary results (with all arithmetic 


carried to 12 significant digits): 


Number of 
intervals Value 
4 250 .52203908 
6 168 .14969914 
8 127 .04780889 
12 86 06241235 
16 6565338565 
pri 45 .35960012 
32 35.29512874 
48 25. 34342609 
64 20 44756141 
96 15 65971496 
128 13.34117182 
192 11.12187172 
256 10 .07933864 


and now our faith in the new tool may be shaken. What 
1s going on? What should be done about it? If the 
next integration we try is not a polynomial, is 1t one 
for which the process works, or is it one like this one? 


These questions invoke predictable answers from 
students: 


1. Take more intervals. 


2. Go to multiple precision--and take more 
intervals. 


If we extend the previous table of results, we find: 


» 
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Number of 
intervals 


512 8.69555371 
700 8. 39606266 
900 8.24235048 
1100 8.15922800 
1500 8.07883197 
2000 8 .03854200 
3000 8.01258448 
5000 8,.00255345 
10000 8.00022129 


We will not do significantly better with multiple precision, 
unless we go to ridiculous extremes (and we have already 
consumed significant amounts of computer time). 


The trouble is due to the nature of the (cleverly 
contrived) function and limits--the function is asymptotic 
to the y-axis. We might do much better if we were to 
break up the integral into two parts: 


from Bee to .1 and from .1 to 1 
Calling these two integrals A and B, we have the 


following results: 


what the problem is, and the weaknesses of 
Never try to substitute brute 


procedure. 
brains. 


With the possible exception of a 


Number of A B sum 
intervals value value value 
4 26 .93022039 2.40790097 29 . 33812136 
6 19 .07052729 2.34178762 21.41231491 
10 12 .95445588 2.31197801 15.2664 3389 
16 9.67450751 2. 30469044 11.97919795 
24, 797010339 2. 30309797 10 .27320135 
32 7.18125527 2.30276351 948401877 
64 6.16264271 2. 30259756 8 46524027 
128 5 .81053812 2.30258590 8.11312402 
256 5.71705426 2.30258514 8 .01963940 
512 5 .69975126 2. 30258510 8 .00233635 
1024 569761527 2. 30258509 8 00020036 
2048 5 .69742899 2. 30258509 8 .00001408 
4096 5.69741581 2. 30258509 8 .00000090 
8192 5 69741494 2. 30258510 8 .00000005 
The main point to all of this is: no numerical 
procedure should be applied blindly. You must know 


the proposed 


force for 
square root 


subroutine, a pathological case can be found for every 
numerical procedure for which it will not work. And 
Elmer's Law says that if you do use some procedure 
blindly, then the next case you try will certainly be 
that pathological case. 


CO 


Problem Solution 


Problem 65 (issue 19) was as follows: 


The 24 odd primes less than 100 are to be placed on 
the 24 faces of four cubes, in such a way that 


(1) 
(2) 


Any toss of the four dice produces a sum 
that is divisible by 4; or 
Any toss of the four dice produces a sum 
that is not divisible by 4 


Are either of these arrangements possible? If the 
24 odd primes are taken to be those from 5 through 101, is 
either arrangement possible? 


If the primes from 5 through 101 are placed on the 


four dice in this way: 


5 53 

13 61 

& Eames 
29 89 

37 97 

41 101 

all of the 

form 4K+1 


then any toss of the four 
and two of the other: 


tf 47 

11 59 

ae ei In any order 
31 79 on each die. 
h3 83 

all of the 

form 4K+3 


Gice will show two of one form 


4K, + 1 
4K, + 1 
AK, + 3 
AK), + 3 


4K + 8 


and the sum is always divisible by 4. 


If the primes from 3 through 97 are used, there will 
be 11 of the form 4K+1 and 13 of the form 4K+3. Then no 
arrangement is possible for either case. (il 
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CONTEST 4 RESULTS 


Contest number 4, 
Square Spiral, appeared on 
the cover of issue 35. 

Consecutive integers 
were to be entered into the 
pattern, beginning as shown 
here, with a square skipped 
after every prime, and two 
squares skipped when two 
primes fell side by side. 
The pattern was to be 
extended and a list of the 
numbers extending north- 
east from the center was 
sought. 


The longest such list 
was produced by Jeffrey 
Shallit, Princeton, New 
Jersey, and is reproduced 
on the following page; 
it has 370 entries, 
indicating that along the way Mr. Shallit kept track of 
over half a million numbers, which included some 45,000 
primes. The zero entries indicate a square along the 
northeast diagonal that remains empty. 


The computing problem involved in this contest is 
not, of itself, very practical. It would be straight- 
forward to attack it using large amounts of storage 
and CPU time. But it could be done, to the limits that 
Mr. Shallit pushed it, with not over 1600 words of 
storage (each at least 20 bits long) and, with careful 
coding, a machine run of perhaps 10 minutes on a modern 
machine. 


In any event, the problem is now useful as a 
coding exercise with known results. 


O 


31786 
32474 
33166 

6) 


34378 
35288 
36013 
36743 
37482 
38232 
38983 
39739 
40511 
41285 
42065 
42847 


0 
4UN56 
45272 
46098 
46924 
47763 
48605 
4ghhg 
50302 
51172 
52036 

0) 


53806 
54701 
55612 
56518 
57422 
58345 


) 
60214 
61165 
62121 
63085 
64049 
65036 
66015 
67006 
67994 
69009 
70030 
71044 


72070 
73099 
74150 
75199 
76251 
77319 
78387 
79468 
80559 
81044 
8272 
83853 
84974 
86103 
87230 
88376 
89526 
90677 
91839 
93003 
94180 
95354 
96544 
97746 


100169 
101390 
102607 
103853 
105093 
106345 
107604 
108866 
110132 
111420 
112709 
114006 
115314 
116623 


0) 
119264 
120601 
121939 
123286 
124639 
125997 
127372 


128748 
130134 
131528 
132928 
134337 
135750 
137168 
138592 
140033 
141473 
142922 
144390 
145857 
147330 
148813 
150291 

re) 


0) 
154796 
156316 
157842 
159387 
160918 
162473 
164023 
165591 
167163 
168753 
170332 

) 


173518 
175125 
176735 
178366 
179979 
181631 
183270 
184919 
18657 

18824 

189911 
191598 


fe) 
194981 
(@) 


198388 
200110 


201837 
203578 
205326 
207070 
208813 
210562 
212336 
214116 

0) 
217703 
219505 
221305 
223128 
224954 
226786 
228622 
230454 
2323i2 

fe) 
236051 
237932 
239813 
241699 
243596 

fe) 


247411 
249327 
251257 
253194 
255135 
257095 
259047 
261023 
26 3005 
264984. 
266968 
268948 
270948 
272958 
274976 
276988 
aiceee 
281068 
283110 
285159 
287233 
289309 


291381 
293473 
295572 
297674 
299792 
301901 
304022 
306144 
308286 
310435 

0 


314736 
316901 
319069 
321241 
323427 
325622 
327821 
330037 
332250 

0 


336719 
338946 
341202 
343453 


0) 
394952 


397 384 
399816 
402275 
4047 32 


(0) 
409666 
412128 
414621 
417111 
419599 
422106 
424606 
427131 
4.29661 
432185 


497957 
500687 
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SR-52 Notes 


The Owner's Manual for the SR-52 programmable 
calculator says (page 75): 


The hait command entered from the keyboard when 
the SR-52 is in the run mode will stop execution 
of a program and return control to the keyboard, 
The program counter is left wherever it happened 
to be at the time of program interruption, 
Program execution will be resumed at that 
location when RUN is pressed, 


This is literally true, but highly misleading. 
It could be taken to mean that while the machine is 
executing a stored program it can be interrupted and 


then restarted at the point of interruption. This is 


not so. The HALT key does not act at the end of a 


command, but at the end of a program step. Thus, such 


simple commands as 


(store the display at register 07) if HALTed at the 
point Indicated by the arrow will restart, if at all, 
with disastrous results. 


Why would one wish to use HALT while a program is 
executing? If the execution of a program takes a long 
time, it might be interesting to monitor it. Or, when 
debugging and testing a complicated program, one might 
wish to run it for a ways and then cut in to see how it 
is doing with the various variables. Or, a program 


known to work properly and give results every two minutes 
‘now has run 20 minutes without results--an interrupt to 


examine storage contents could reveal troubles (which 
might not exist--the data has changed) and a RESTART 
could salvage the 20 minutes of calculation. 


There can be many legitimate reasons for using a 
HALT on a running program. To be sure, if the machine 
has the printer attachment, most of these situations can 
be covered by suitable PRINT commands inserted within 
the program, but the printer sells for almost as much as 
the calculator itself (currently $250 for the printer 
vs. $300 for the SR-52). The manual speaks in glowing 
terms of the virtues of having a printer. 


There is a technique for inserting a sense-switch 
HALT in a program, provided that the program uses no 
trig functions. The SR-52 is set to radian or degree 
mode by means of an external slide switch. If the 
switch 1s set to radian mode, then the calculation of 


(SINE 30 - .5) 


will be non-zero. This can be tested in the program, 
and a HALT can be conditionally programmed which will 
take effect during execution only when the switch is 
set to degrees. 


We seem to have neglected the consecutive numbering 
of Problems as they have appeared. The following list 
will bring the system up to date: 


Problem 
number Name Reference 
126 Life or Death PC37-14 
127 Peripatetic Jumping Bean PC38-2 
128 K-level Sieve PC38-7 
129 Circuitous Race PC38-18 
130 Ring-a-ding PC39-1 
131 Outguess PC39-4 
132 8 Dice PC39-17 C] 
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Summing 7-Card Decks @ 


A deck of cards bears six 3-digit numbers on each @ 
card, in columns 1 through 18. It is formed of a number 
of seven-card sub-decks; each sub-deck is identified by 
the 3-digit number in columns 1-3 of the first of the 
seven cards. 


As the cards are read, a total is to be formed of 


the other 41 numbers in each sub-deck. After each 9 
sub-deck has been handled, its identifying number and td 
the sum for that deck is to be printed. oa 
= 
Assume that a subroutine is available that will a 
read a card and place its six numbers in words addressed 2 
at G, G1, G+2, G+3, G4, and G+5. The ldentifying 
numbers for the sub-decks are all over 500; all data 
numbers are less than 500. The end of the full deck 
is signalled by a sentinel card bearing the number 999 
in its first three columns. 
The following two sets of seven cards would produce 
the printed lines shown at the bottom: 
‘557 001 002 003 OO4 005 
006 007 008 009 O10 O11 
012 013 O14 O15 O16 O17 
oF 019 eee 021 Bee 023 
peo oe ee 0) 863 100 101 101 102 102 
030 031 032 033 034 035 103 103 104 104 105 105 
036 037 038 039 O40 O41 406 106 107 107 108 108 


109 109 110 110 i111 iit 
112 #112 #113 #113 «114 114 
11 11 nae; alae; ayy ay 
11 11 119 119 120 120 


557 861 
863 4520 


A flowchart for a possible solution to the problem 
as stated is shown. 


This problem is of more than passing interest--it 
will be referred to in later issues. Notice, in the 
logic of the proposed solution, that provision has been 
made for the case in which a sub-deck of less than 7 
cards appears. 


LT-On0d 


(oa 

(S+0)+( 4+) 
+(E+0)+(2+0)+ 
(T+0)+(0)+(S) 


(s) “(ano ) 
qutad 


(s) ‘(in0) 
‘aZessow iota 


quTad 


*spareo syuno0d 9g aaqgunog 
*S UT paWwsoy ST ums partnbay 


“G+D SHt+D fE+D 

‘Z+D ‘T+D “D SpuoM ayy OF 
pPaeo suo Jo sguaguod 3u4 
sgnduy aqyayw euTyZnorgns au, 


swoad auvo-2 ONIWWAS 


(Z| (s) — (s+0)+ 


(7+0)+(€+0)+ 
(2+))+(T+9)+(S) 
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Problem Solution 


Te Stack Action Problem (PC37-3) called for a 
subroutine to stack numbers (which are limited to the 
range 001 to 999) in ascending order as they arrive from 
the main routine. No number is to be put in the stack 
if it lies within 10 of any previously stacked number, 
and the procedure ends when 50 numbers have been stacked. 


The APL code given here comes from M. E. Sandfelder 
and @. Truman Hunter, Poughkeepsie, New York. The 
pase part of the code is found on lines 2, 8, 12, 
and 14. 


The problem statement in issue 37 called for a 
test procedure, which is still needed. 


VY R*BUILD X 
[1] aINITIALIZES RESULT R TO THE VALUE 1000 
(2) R+1000 
[3] ad ONE LINE LOOP THAT PICKS A RANDOM NUMBER UP TO 999 ASSIGNS IT TO 
C4] aA VARIABLE N SUBTRACTS IT FROM ALL VALUES OF R TAKES THE MAGNITUDE 
[5] aOF THE RESULT TESTS TO SEE IF ALL DIFFERENCES ARE GREATER THEN 10 AND 
[6] aIF SO BRANCHES TO NEXT LINE, OTHERWISE STAYS ON LINE 8 AND REPEATS 
(7] AOPERATION. 
- [8] +(V/102|R-N+?999)/8 
[9] aACATENATES N TO THE RESULT R TESTS TO SEE IF THE NUMBER OF ELEMENTS 
(10] alN R IS LESS THEN THE REQUIRED NUMBER CONTAINED IN THE DUMMY VARIABLE X, 
[11] aAND IF SO GOES BACK TO LINE 8 ,OTHERWISE BRANCHES TO NEXT LINE OF PROGRAM 
[12] +(X>pR«R,N)/8 
[13] ARBARRANGE RESULT R IN ASCENDING ORDER USING GRADEUP AND INDEXING 
C14] R+RCAR] 
Vv 
15] V 


RR*+BUILDN 50 
5 10pRR 


3 34 71 91 104 131 142 154% 4179 £4200 
213 233 252 270 287 307 318 351 365 386 
412 427 439 455 468 498 513 537 555 572 
591 625 652 681 702 714% %725 748 %783 802 
817 846 860 876 888 901 937 9355 983 1000 
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Proceed from the upper left cell 


to the lower right cell (marked B) in 


lines, with each le 
have the sum of the contents of the cells traversed 


over 6 cells long, with the path not crossing and 
the greatest possible. 


In the pattern shown, there are 256 cells, each 
not touching itself anywhere. 


containing a 3-digit number at random. 
such a pattern was the basis for news 


under rules such as these: 


The Gemeroy Problem 
ago, 


A possible path is indicated in the 
figure, with a total of 69100. 

It should be possible to find paths 
with larger totals. Notice that 
the rule about the path not touching 
itself rules out paths like this: 


PC40-20 


Finding a path with a larger total may be easy. 
The difficult problem is a systematic attack to find the 
path with the largest total. The rules for getting from 
A to B are sufficiently stringent that the number of 
possible paths is not very large. Given the 256 cell 
values in storage, then each path that can be developed 
geometrically can be applied four different ways to the 
array, and the required sums can be obtained. The 
trick 1s to arrange for the controlling program to make 
the minor alterations on the paths. For example, in 
the path shown on the array, from the cell containing 
the value 867 (row 11, column 16) there are at least 6 
other paths possible to point B. 


The newspaper contests have disappeared, probably 
because people started using computers on such problems. [e] 


NOQNONNDONDONO NNO NNO NDNNOONONON0000 000 
H one programmer can daa tack in one day, 


two programmers can da it ta two days 
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